1 Introduction

The permafrost landscape of the Arctic shoreline is one of the most dramatically changing environments in the world. Despite the fact that Arctic permafrost coastlines represent 34% of Earth’s coast and are affected by permafrost [1], the physical processes and mechanisms that drive change along these coastlines remain poorly understood [2, 3]. Along the Arctic shoreline, coastal erosion is expected to dramatically increase over the next decades due to the combined effect of longer and warmer thawing seasons, declining sea-ice cover that allows larger waves and storms to collide with the coastline at both longer durations and at higher coastal elevations. With this heightened erosion one can expect that constituent fluxes both into the sea and air will also increase over time.

Erosion rates appear to be increasing along some sections of the Arctic coast [4, 5] reaching up to 25 m/year [6]. This has resulted in an annual efflux of 14 Tg of particulate organic carbon into the coastal ocean [7]. This carbon flux is of equal magnitude as the annual delivery from all Arctic rivers, or the vertical net methane (CH4) emissions from terrestrial permafrost [8, 9].

While it is known that thawing permafrost creates up new areas for methane emissions, the identification of these pathways using traditional ground-based or satellite-based methane-detecting methods has been both challenging and insufficient [10]. In coastal permafrost environments, pathways from which methane can be released into the atmosphere [11] are typically (a) the top half meter of soil that thaws each summer (active layer), especially from thermoerosional drainages or meltwater lakes where microbes are most active; (b) escape conduits from disintegrating permafrost caps and/or ice wedges; and (c) from thawing offshore permafrost and gas hydrate dissociation (Fig. 1). Here we focus our efforts on the first of the above three pathways, where high methane fluxes are the result of thawing permafrost that is associated with saturated systems such as thermokarst bogs, fens [12,13,14] and as in this study thermoerosional gullies [15, 16].

Fig. 1
figure 1

Three common pathways for methane release in coastal permafrost bluff environments

Spatial identification and quantification of CH4 efflux hotspots may serve a twofold purpose along coastal bluffs. First, quantifying CH4 efflux from hotspots is essential for accurate methane budgets. Secondly, these hotspots could also indicate areas of lower soil stability, thereby identifying coastal sites highly vulnerable to erosion. Identification of methane pathways may therefore help gain better local insights on permafrost integrity and help with erosion prediction models.

Currently, spatial methane distribution in the Arctic is being recorded mostly with top-down methods by several different means—all of which come with major drawbacks: (1) tower station measurements have a sparse spatial coverage and cannot identify individual methane sources; (2) satellite measurements can reveal large scale atmospheric background signals but give little indication on the ground-based sources of methane sources; (3) airplane-based measurements can only be carried out at a minimum height of ~ 100 m and are not as accurate as in situ data. All of the above methodologies are associated with significant costs unattainable for most institutions and stakeholders (e.g. [17]).

Until recently, methane measurements with systems sensitive to parts-per-billion levels were also confined to large, heavy and expensive laboratory equipment; now, a new class of significantly more cost-effective and lighter tunable diode-laser using absorption spectroscopy [18] allows for the measurement of methane using small Unmanned Aerial Systems (UAS). While few of these methane-detection systems have been developed (i.e. [19, 20]) none have aimed to create an off-the-shelf, cost effective UAS design. An overview of the challenges and limitations of using small UAS systems to monitor methane and other atmospheric gases has been compiled by Villa et al. [21]. While this manuscript focuses on the application of an unmanned methane sensor in a very remote, largely undisturbed location north of the Arctic circle, where permafrost is thawing, the proposed system could also be utilized for more urban environments and for pipeline leak detection. Emran et al. [22] for example, have recently demonstrated the feasibility and advantages of UAS-based methane monitoring for grass-topped landfill degassing, and Bretschneider et al. [23] showed how such a system could detect leaks in pipelines. The methane laser sensor used herein was first developed in 1992 [24] and became commercially available in 2013 [25]. The sensor has been shown to function in a natural environmental setting to demonstrate time variation of low CH4 concentrations over rice paddies [26] or small ponds and ditches [27]. We here present the design and components of a CH4—UAS system.

In this paper, we demonstrate our efforts to detect methane being emitted from thawing coastal permafrost and map these concentrations. Our study was carried out along the coastal bluff environment of Barter Island, AK, which is part of a larger Arctic coastal erosion investigation study by the U.S. Geological Survey (USGS; i.e. [28,29,30]).

2 Methods

2.1 The UVA-platform

In this study the 3DR Solo was used as the Unmanned Aerial Vehicle (UAV) platform (Fig. 2). It was chosen because it is currently the only small UAV platform that has been approved for use by the USGS for research purposes by the Department of the Interior. Although not ideal due to its low payload of 700 g the rotary-wing quadcopter is inexpensive (< $500 USD) and was upgraded with a highly accurate professional grade GPS system utilizing European (Galileo), Russian (GLONASS) and North American (GPS) positional information. The GPS system has a Horizontal Dilution Of Precision (HDOP) of 0.6 inside the UAV [31]. At full payload and at 3–6 °C, the maximum achieved flight-time was between 8 and 10 min. For future flight configurations, we highly recommend utilization of a UAV platform that has a greater payload and a longer flight time, but similar GPS accuracy and flight stability.

Fig. 2
figure 2

Required UAS components with weights: a Methane sensor (530 g), b micro Android OS (55 g), c UAV (maximum payload 700 g) with UAV leg extensions (60 g), d vibration mount (20 g), e mount for sensor (5 g)

The vibration dampening methane sensor mount was 3D-printed (Fig. 2) and attached to the central Accessory Bay, the area behind the gimbal under the 3DR Solo. It is intended for secondary accessories, including additional communications hardware and other high power devices [32]. Total cost for UAS and 1 month sensor rental was < $2500.

2.2 The methane sensor

The Pergam Methane mini-G (SA3C50A) is currently the only commercially available methane sensor that has low enough weight and high enough range (1–50,000 ppm) and sensitivity (1 ppm) as well as a fast enough recording speed (10 Hz) to allow for remote sensing use (Table 1). The sensor uses Infrared Absorption Spectroscopy and a semiconductor laser for methane detection [18] and is advertised by the manufacturer for aerial use over non-flat surfaces such as grass-covered landfills. While the strongest absorption bands for methane exist at ~ 3300 nm and ~ 7700 nm these bands are overlapped by absorption bands of water vapor and CO2. In order to avoid this overlap, a weaker band frequency of 1653 nm has been determined to be free of overlap by absorption bands of water vapor and CO2 (see “Appendix”, Fig. 8). Nonetheless it should be mentioned that tunable laser absorption spectroscopy is an ongoing research area and others [33] have made clear that overlap of bands can occur in field observations if the sensor is not finely tuned to the specific band. Consequently, the detector relies on transmitting a laser beam with a frequency tuned precisely to the absorption characteristics of methane gas (1653 nm) [34]. The integrated concentration of methane between the sensor and the target point is measured by transmitting a detection laser beam towards the target point, and then detecting a fraction of the diffusely reflected beam from that target point. The measured value is expressed by an inventory (ppm*m) which is the methane concentration (ppm) multiplied by the distance (m) between the sensor and its target. All measurements were normalized to 1 m distance to reflect the more common volumetric ppm values found in atmospheric sciences literature. Consequently, the measurements presented herein represent path-averaged concentrations in ppm, which implies that ground concentrations could feasibly be higher (Fig. 3).

Table 1 Methane sensor technical information
Fig. 3
figure 3

UAV-based tunable diode laser spectrometry. The system calculates path integrated methane concentration in ppm*m. This value is then divided by the distance to the reflection surface to get path-averaged concentration in ppm which represents the more common notation standard (e.g. [35]). In example A, this procedure ((150 + 2+2 + 2+2)/5) yields a concentration of 31.6 ppm, while B ((2 + 30 + 2+2 + 2)/5) has a concentration of 7.6 ppm. In the illustration ~ 2 ppm/m was chosen to present a typical CH4 background concentration

Methane concentration measurements were updated at a rate of 10 Hz. In addition to the internal calibration that the instrument undergoes during each reboot, a calibration experiment was carried out using quartz glass calibration gas cells at varying distances and concentrations. While the manufacturer lists the CH4 levels with a ± 10% accuracy our calibration experiment suggested an accuracy of ± 30–40%. Calibration curves are shown in the “Appendix”. All data were GPS referenced and stored via a standard Bluetooth link on a UAV mounted Android OS device (Fig. 2).

2.3 Flight plans and post processing

Automatic flight paths were created using Mission Planner [36] and executed to fly at a continuous speed of 3 m/s. A high precision Lidar-based Digital Elevation Model (DEM) [15] was used to guide the methane sensor at a constant distance of 5 m above ground. The normalization of the data to 1 m distance occurred by simply dividing the set flight elevation by five. The flight elevation was chosen in part because the sensor is most precise at distances of 5–10 m (see “Appendix”) and because it was the minimum height at which the drone could be safely operated above potential human activity. However, there was no recorded human activity nor engineering construction that could have affected the measurements within the study site. The flight paths were created to be above the flat vegetation tundra landscape, slightly inland (~ 5 m) of the coastal bluff (Fig. 5). All flights were carried out during zero to low (< 3 m/s) wind conditions and flown against the prevailing western wind direction. All human infrastructure that could potentially act as an anthropogenic methane source was downwind from our study site. This strategy was chosen in order to identify individual hotspots along the flight path rather than following a plume. Air temperatures during all flights were at a constant 2 °C and between 35 and 45% humidity as determined by the nearby (located 650 m away) airport weather station. The methane sensor automatically marks data that was recorded with a below minimum reflected energy as containing an error such as over bodies of water. This data was filtered from the utilized dataset following a previously published MATLAB (R2007a) methodology described by Emran et al. [22]. Visualization of the methane data was carried out using Google Earth software but can be performed using any GIS software program.

3 Results and discussion

Field deployment of the UAS was carried out on Barter Island, Alaska between Sept. 3–5, 2017 over a near linear path of 2 km closely following the coastal bluffs. In addition to the previously performed calibration experiment in the lab (see “Appendix”), a test grid was flown at a higher altitude of 10 m over an area where permafrost was expected to be largely intact and methane levels relatively constant (Figs. 4a, 5). The purpose of this test grid was to demonstrate the functionality of the sensor in situ and establish a CH4 atmospheric baseline as well as background noise. Because atmospheric background levels are very close to the sensor’s sensitivity rating of 1 ppm, a higher altitude was chosen to allow for a greater path-integrated measurement in the background level test. The background level test flight revealed a relatively high instrument variance of 1.17 (± 56%) over 5600 measurements but also a consistent linear average of 2.06 ppm (Fig. 4a). Considering the manufacturer’s reported error range (± 10%), this value closely mirrors the average regional background measurements of 1.95 ppm for Sept., 2017 (Fig. 4b) established by the NOAA’s Global Monitoring Division [37].

Fig. 4
figure 4

a (Lower) Methane background levels from Test Grid flight (see Fig. 5 for location) were recorded as ppm × m measurements and normalized to ppm by division of the flight height of 10 m. b (Upper) Long term regional atmospheric background methane signal

Fig. 5
figure 5

North facing bird’s eye view (45-degree) of methane hotspots and major melt-water pathways in a coastal permafrost environment on Barter Island, AK. The 2 km long UAS path starts at 70° 7′56.39″N 143°42′1.28″W and ends at 70° 8′1.05″N 143°38′49.79″W

Our first results from Barter Island, AK show that spatial methane distribution along the 2 km coastal path is mostly at atmospheric background levels of 2 ppm ± 35%. There are however several large methane peaks within the flown path. These peaks range between 20 and 210 ppm and are spread out along the path. After closer analysis of the surface topography it became clear that the peaks were linked with areas of high permafrost thaw and thermoerosional gullies (Fig. 5).

High methane fluxes from thawing permafrost are often associated with saturated systems such as thermokarst bogs, fens, and lakes or low centered polygonal tundra [12]. Much less research has been focused on thermoerosional gullies because they may not typically be saturated and they may represent a small fraction of the landscape. However, recent studies [15] have shown high CH4 fluxes from gullies such as ours, even if soil moisture is not changed relative to the undegraded landscape. The mechanisms driving these fluxes and the contribution of these fluxes to the overall coastal methane budget, are further issues to be resolved.

A comparison between methane hotspots and long-term coastal bluff retreat (Fig. 6) showed that the highest methane peaks could also be associated with some of the highest retreat rates along the coastal bluff. More specifically, the methane hotspots mark thermoerosional gullies, which are also areas of amplified permafrost thaw and thus can be seen as erosional hotspots along the bluff coastline. However, not all peaks correlate, which indicates that while thermoerosional gullies can be areas of lower bluff stability, the primary bluff retreat is likely caused by a combination of thermo-denudation and thermo-abrasion [38].

Fig. 6
figure 6

Comparison between coastal bluff retreat rates (blue line) [6, 38] and methane concentrations (red line)

The main uncertainties associated with this study are caused by potential variations in the ground to sensor distance. These variations are most likely associated with non-linear movements of the UAS such as roll, tilt and yaw that are produced by flight turns and variable wind speeds and directions. While the data on off-nadir angle could not be recovered due to an overwrite, common deviations from Nadir for this UAS-class are < 15° [39, 40]. Assuming a 15° deviation from Nadir, this would lead to a maximum 3.5% increase in measured methane levels, which is significantly below the published accuracy (± 10%) or our measured accuracy (± 35%) of the instrument (Table 1, “Appendix”). Nonetheless, in order to account for this effect we included many turns in our test grid at constant flight speed (Fig. 4a) forcing yaw, tilt and roll and thus establishing a conservative measure of instrument background noise. Conversely, in the permafrost bluff survey we attempt to minimize this effect by avoiding all turns in our flight path (Fig. 5). Unfortunately, wind, fog and snow prevented additional flights during the time-limited field research campaign, and locations north of the Arctic Circle such as Barter Island are not easily accessed to allow for low-cost follow up studies. Nonetheless, with respect to future deployments, we hope to employ traditional ground-based point-sampling techniques to better quantify the relative limitation of each methodology in more detail. Additional uncertainties in our measurements may be the result of the under or over represented gullies in the DEM used in this study. Gullies ranged in depth by 1–5 m and an additional conservative ± 20% measurement error can be estimated from comparisons between field elevation observations of gullies and DEM results. While this leads to a total uncertainty of ± 55%, which is mirrored by our field observed uncertainty of ± 56%, the reported peaks found in gullies are more than 2000% above background levels (Fig. 5) and thus remain significant with these uncertainties applied.

Despite the current uncertainties, the herein presented UAV-based methodology falls inside an important spatial gap between local and regional scales that was previously difficult to bridge by existing methodologies. Methane measurements occur along a range of spatial and temporal scales (Fig. 7), from small-scale near-instantaneous measurements of emissions from individual sources (such as from a smokestack) to large-scale global assessments of annual emissions (such as from satellites). Associated with the different methane sensing technologies are major strengths and weaknesses that contribute to the applicability and scope of the produced data (Table 2).

Fig. 7
figure 7

Spatial and temporal scales associated with major methane sensing platforms and methodologies (e.g. [17])

Table 2 Advantages and disadvantages of major techniques measuring methane emissions

The herein presented technology could also be used to look for areas of geological methane release and CH4 release from methane hydrates. However, in those cases, isotopic methane composition will help in answering questions of geologic methane origin. Additionally, planned future studies focusing on the onset of the thawing season are likely to answer the question if the herein presented UAS methane-detection methodology is capable at predicting erosion hotspots through early identification of methane hotspots.

Spatial methane distribution data that can be used to identify source pathways is currently still very rare and traditionally its monitoring has been reserved to large institutions capable of absorbing the associated high costs. Compared to previous efforts that collected spatial methane data in the Arctic by helicopter [41, 42], by small aircraft [10], by ground based laser scanning [43], or with a ground based solid-state trace gas sensor [44] the methodology presented herein exemplifies a significantly more affordable (< $2500) and more flexible method towards acquiring spatial methane distributions.

With the rapid commercial spread of UAVs the presented UAS system has the potential to accelerate the availability of high resolution spatial methane data through university and citizens science and thus gain significantly better insights into methane pathways and its related processes.